Breakdown of Kolmogorov scaling in models of cluster aggregation with deposition 
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The steady state of the model of cluster aggregation with deposition is characterized by a constant 
flux of mass directed from small masses towards large masses. It can therefore be studied using 
phenomenological theories of turbulence, such as Kolmogorov's 1941 theory. On the other hand, the 
large scale behavior of the aggregation model in dimensions lower than or equal to two is governed 
by a perturbative fixed point of the renormalization group flow, which enables an analytic study 
of the scaling properties of correlation functions in the steady state. In this paper, we show that 
the correlation functions have multifractal scaling, which violates linear Kolmogorov scaling. The 
analytical results are verified by Monte Carlo simulations. 

PACS numbers: 05.10.Cc,47.27.Gs,05.70.Ln,68.43. Jk 



Understanding Navier-Stokes turbulence is a big chal- 
lenge of modern theoretical physics. It is a nonlinear 
system far from equilibrium with no obvious small pa- 
rameters that would allow a perturbative treatment. It 
would therefore be instructive to study simpler nonequi- 
librium models which possess the qualitative features of 
turbulent systems, yet are analytically tractable. Such 
an approach has proved to be fruitful, especially during 
the past decade. The study of Burgers turbulence 0, the 
Kraichnan model of passive scalar advection and turbu- 
lence in kinematic magneto hydrodynamics 0] led to a 
better understanding of the limitations of phenomeno- 
logical theories of turbulence. Also, it concentrated at- 
tention on more general concepts of turbulence, such as 
intermittency and breakdown of self-similarity. 

In this paper we study breakdown of self-similarity in a 
model of diffusing, coagulating masses in the presence of 
a steady influx of monomers. Consider a ci-dimensional 
hypercubic lattice whose sites are occupied by particles 
that carry positive masses. Multiple occupancy of a site 
is allowed. Given a configuration of particles, the system 
evolves in time as follows. At rate D, a particle hops 
to one of its 2d nearest neighbor sites. At rate 2A, two 
particles on the same site coalesce together to form a 
new particle whose mass is the sum of masses of its con- 
stituents. At rate J/ mo, a particle with mass too is in- 
jected at a site. The parameter J is the average mass flux 
into the system. For simplicity of analysis it is assumed 
that the rates D and A do not depend on the particles 
masses. The initial condition is one in which there are 
no masses. We call this model the mass model (MM). A 
generalization of this model is one where the masses (or 
'charges') could be positive or negative. Here, with rate 
Jc/toq particles of charge too and —mo are input into the 
system. J c has the meaning of average influx of square of 
the charge. This model will be called the charge model 
(CM). We will be interested in the continuous limit of 
these models. 



A feature of the steady state of MM is the presence 
of constant flux of mass from small masses to large 
masses via coagulation. This is analogous to turbulent 
systems where there is a constant flux of energy from 
small wave numbers to large wave numbers via nonlin- 
ear interactions. Let C n (m\, . . . , m n )(AV) n ]X drrii be 
the probability of having particles of masses in the in- 
tervals [mi,m,i + drrii] in a volume AV^. We ask how 
C„(mi, . . . , m n ) varies with mass when mi, . . . , m n 3> 
too- In particular what is the value of the homogene- 
ity exponent j(n) defined through C„(Ftoi, . . . , Tto^) = 
r~ 7 (")C n (TOi, . . . , to„)? The quantity in turbulence that 
is analogous to C„(toi, . . . , to„) would be (E^ ■ ■ ■ E^ n ), 
where E^ is the energy corresponding to the wavenum- 
ber hi. 

MM has been studied in many different contexts. Ex- 
amples include submonolayer epitaxial thin film growth 
by deposition of atoms onto a substrate in the limit when 
the distance between clusters is much larger than typical 
size of a cluster jjjj, river networks 0, Kg, force fluctu- 
ation in in granular bead packs and nonequilibrium 
phase transitions 0, H] ■ It was one of the first models for 
self organized criticality wherein power laws are gener- 
ated from simple dynamical rules. In this context, it also 
maps 9] onto the abelian directed sandpile model of self 
organized criticality 01 ■ Finally, the mean field limit of 
MM is mathematically similar to the kinetic equations of 
three- wave weak turbulence . 

It was shown in Refs. [l2L Hij . fl4| that in one dimension 
7(1) = 4/3 for MM and 7(1) = 5/3 for CM. By studying 
the two point correlations, it was shown that in d < 2, 
7(1) = (2d+2]/ (d+2) for MM and 7(1) = (3d+2)/(d+2) 
for CM 

In this paper, using the renormalization group (RG) 
formalism, the exponent 7(71) is calculated as an expan- 
sion in e = 2 — d up to order e. In two dimensions, the 
upper critical dimension of the model, the logarithmic 
corrections to the mean field results for C n {m\, . . . , to„) 
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are calculated. Exact results are obtained for C^(mi, 7712 ) 
in all dimensions. The Kolmogorov prediction for 7(71) 
based on self similarity is shown to break down. 

We first determine the dependence of C n (m\, . . . , m n ) 
on mass using a self-similarity conjecture similar to Kol- 
mogorov's 1941 conjecture about the statistics of velocity 
increments in hydrodynamic turbulence. Assume that 
C n depends only on the masses rrii, mass flux J and the 
diffusion coefficient D. The dimensions of the various 
parameters describing the continuous limit of the model 
are [J] = ML~ d T~ l , [D] = L 2 T~ X , [C n ] = L- nd M~ n 
and [to] = M. There is a unique combination of D, 
to and J that has the dimension of C n given by C n ~ 
(JI}- 1 )«<V(<*+ 2 ) m -'rfcoim(n) ; w here 



lkolrn{n) 



2d + 2 
d+2 



n. 



(1) 



is the Kolmogorov scaling exponent. As expected, the 
dependence of jkoim on the index n is linear, reflecting 
the assumed self-similarity of the statistics of the local 
mass distribution N(x,m,t). When n = 1, ^koim 
' + 2) /(d + 2), which agrees with the result of an exact 
computation for d < 2 [16J. 

The self-similarity conjecture assumes that C„ does 
not depend on the following: the reaction rate A, the 
lattice spacing, the position of the source m and the 
box size AVdmi . . . dm n . The lack of dependence on the 
lattice spacing is expected due to the renormalizability 
of the effective field theory describing MM below two 
dimensions. We will however find an anomalous depen- 
dence of correlation functions on a length scale depend- 
ing on the other parameters that leads to the violation 
of self-similarity. 

Starting from the lattice model, it is possible to con- 
struct an effective field th eory of MM using the formalism 
due to Doi and Zeldovich |l7lll8j. Furthermore, it is pos- 
sible to establish an exact map between this field theory 
and the following stochastic integro-differential equation, 

MM. 



8 \ r m 

^- - L>V 2 cj)(m) = A / dm'(j)(m')(j)(m - to') 
9t J Jo 

- 2\(f)(m)N + —5(m - to ) + iV2X(j)(m)r](x , t)(2) 
toq 



where N = f Q dm' cf>(m') , i = \J— 1, and i](x,t) is 
white noise in space and time with (r](x, t)r](x' , t')) — 
5(t — t')5 d (x — x 1 ). All correlation functions of the mass 
distribution can be expressed in terms of the correlation 
functions of </>(m, x, t) In particular, 



C„(m) 



1 



(m,x, t)} 7 



(3) 



where (...) denotes averaging with respect to noise rj \2l\ . 

Without the noise term, Eq. reduces to the mean 
field Smoluchowski equation of the model. Thus, all the 
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FIG. 1: Propagators and vertices of the theory. 

fluctuation effects are encoded in the imaginary multi- 
plicative noise term. Equation @ simplifies after takin 
Laplace transform with respect to the mass variable [2 
Let R^(x,t) = N(x,t) - J Q co dm(f)(x,m,t)e-^ m . Then, 

'd 



at 



\Ri + — + W2XR^(x, t), (4) 



too 



where j — J(l — e^™ 10 ), and the dependence of on 
(x, t) has been suppressed. 

Equation Q has the form of the stochastic rate 
equation of the A + A — > A reaction in the pres- 
ence of a source [2^], reducing the the computation of 
the mean mass distribution in MM to solving a one- 
species particle problem. However, in order to com- 
pute C n (m,t), the correlation functions of the form 
(x, t)R^ 2 (x, t) . . . R^ n (x, t)) need to be determined. 
These are non-trivial, as the stochastic fields R^(x, t)'s 
are correlated for different values of \x via the common 
noise term in Eq. (|4"|). 

The set of Feynman rules for perturbative computation 
of correlation functions C n follows from Eq. Q |23j and is 
summarized in Fig. ^ The n-point correlation function 
(niLi R-Hi (%ii i*)) is given by the sum of all Feynman 
diagrams that have n outgoing lines built out of blocks 
shown in Fig. ^ 

Let R m f, denoted by a thick line with a cross, be the 
sum of all tree diagrams with one outgoing line. The 
equation satisfied by R m f is shown in diagrammatic form 
in Fig. dR m f/dt = j/m — \R 2 m j- This corresponds 
to the noiseless limit of Eq. Q. Solving, 



R m f(t) 



J 



mo A 



tanh 




(5) 



To account for the noise term, we have to include di- 
agrams with loops. To construct the loop expansion, 
it is convenient to introduce a tree level propagator 
G m f(x2t2',xiti). The equation obeyed by it is shown 
in Fig. |33. The solution is 



Gm/(2;1) 
G (2;l) 



cosh 



cosh - 



-fi(t 2 -ti) 



(6) 



(A) 



-X = 



(B) 




FIG. 2: Diagrammatic form of mean field equations for: (A) 
Rmf and (B) i?_R-response function. 



where Go is the Green's function of the linear diffusion 
equation, 2 = (£2, £2), 1 = (£i,ti), and fl = 2y/ j\/mo 
is the inverse of the mean field response time. All terms 
in the loop expansion constructed using the vertices of 
FigH G m f and R m f are finite in d < 2. 

We now show that the loop expansion corresponds to 
weak coupling expansion with respect to A. Consider a 
diagram contributing to (R 1 ^) that has L loops, V vertices 
and N i? m /-lines. The A factors arise from Ld momen- 
tum integrals (A id / 4 ), V time integrals (\~ v/2 ), N R mf 
lines (\~ N / 2 ) and V vertices (\ v ). Thus, the correspond- 
ing Feynman integral is proportional to \- N / 2 + v / 2 + Ld / A . 
Also, the number of triangular vertices in the graph 
equals N — n and the number of quartic vertices is equal 
to the number of loops L, i.e., V = L + N — n. Thus, any 
i-loop graph contributing to the mean mass distribution 
is proportional to \~i+ L (H i~\ We conclude that loop 
expansion corresponds to the perturbative expansion of 
(i2™) around i?™^ with the parameter A^^. 

The conditions under which the loop corrections to 
the tree level answer can be neglected may be derived 
using dimensional analysis. The scale of diffusive fluctu- 
ations is given by the only constant of dimension length 
which can be constructed out of D and and j /m®: Lp = 
(mo-D/ ^j) 1 /( d + 2 ) . The dimensionless expansion parameter 
in the loop expansion above is therefore g{Ljj) = XL e Dl 
where e = 2 — d. The large mass behavior of C n (m) is 
determined by the small-/it behavior of (i?™). In d < 2, 
go — ► 00 when \l — > and the loop expansion breaks 
down. Thus, a re-summation of loop expansion is needed 



to extract the fj, — ► behavior of R^. This will be done 
using dynamical RG formalism 




FIG. 3: Zero and one loop Feynman diagrams contributing 
to {R l _ ll (x,t)R fl2 (x,t)) . The circles mean that all outgoing 
lines terminate at the same spatial point. 



depend on the length scale Lq. Under RG, gn goes to 



Large Ld 



g* as L D — > 00 [23I . Hence, 
corresponds to small /i. Then j/mo ~ J/i, and 
(Jfj,) d /( d+2 \ The inverse Laplace transform gives 



Ci(m) 



m-i jW/(«j+2) 



d < 2, 



(8) 



TO 7(1) 

where 7(1) = (2d + 2)/(d + 2). This result is exact to all 
orders in e. 

The RG analysis shows that C\ depends only on J, m 
and D and thus verifies the self similarity conjecture for 
n = 1. Therefore, it is not surprising that 7(1) coincides 
with 7fc j m (l). 

We now examine higher order correlation functions of 
the form {R tll {x 1 t)R^ l2 {x,t) . . .). Now, extra ultraviolet 
divergences, that do not get canceled by coupling con- 
stant renormalization, appear in the perturbative expan- 
sion of these composite operators |25j. These extra di- 
vergences lead to a spectrum of anomalous dimensions 
which are calculated below. 

The diagrams contributing to {R^iR^) up to one loop 
are shown in Fig. [3] It is straightforward to generalize 
them to n > 2. For n > 2, diagrams contain n outgoing 
lines and Q) connected one loop diagrams. Computing 
the diagrams using equations J2J) , JUJ and simplifying, we 
obtain 



We first examine (R^). It was shown in Ref. [13, E^] (TT-R/ii) =TT(-Rp 



that renormalization of the coupling constant alone reg- 
ularizes the loop expansion for when e — > 0. In par- 
ticular, (Rfj) has no anomalous dimension and therefore 
has the form Ljf "f[gn,Lo/ 'Ljj] where Lq is a reference 
scale, gn is the renormalized reaction rate and is related 
to g(L ) by 



i=i 



i=l 



n(n - l)Ar(e/2) 

(87TL>)' i / 2 (f7l+f72) e / 2 " 



(9) 



Expressed in terms of the renormalized coupling gn, 
Eq. I© reduces to 



(n^)=n^(^)) 



9R = 



g(Lo 



(10) 



l+g(L )/g*> 
where g* = (87r) d / 2 [2r(e/2)]- 1 [Uli^. 



(7) 



does not 



Since coupling constant renormalization removes all sin- 
gularities in (.R/j (<?#)), there are no singularities in the 
product n"=i {Rfx,i{9R)) when e — > 0. 
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The remaining singularity is canceled by multiplica- 
tive renormalization of the composite operators. Let 
(11™= l r m)r = z n(W=i R^i), where Z n is fixed by the 
condition that (11™= 1 R^i) R 1S n °t singular in e 26]. Thus 
Z n = l + n(n-l)g R /{4ire) + 0(e 2 ). The bare (fl^i 
is independent of the reference scale Lq. Thus, 



d 



Z n 1 (W R ^)r 



= 0. 



(11) 



Let Li — [moI?/j(/Xi)] 1 ^ £i+2 - ) . By dimensional analysis, 
(Yli=i R ^) = L o 2d f( L i/ L o,9R), where / is a dimen- 
sionless scaling function. Substituting into Eq. (|llfl. and 
solving the resulting equation in the limit Li 3> Lq, we 



obtain (fliLi Rpi 



is a scaling function of the variables Li/ Lj. 
Laplace transform gives 



, where $i 
An inverse 



C„(mi, . . . ,m n ) 



n 1 

$ 2 n— 



J 



mi \Drrii 



rf + ,:(n-l)/2 
d + 2 



(12) 



1UW 

one loop 



FIG. 4: The variation of homogeneity exponent 7(71) with n 
is shown for one dimension. The dotted line corresponds to 
Eq. 11311 with e = 1 and terms of order e 2 and higher set to 
zero. The values of 7(0), 7(1) and 7(2) are exact while 7(3) 
and 7(4) were obtained by Monte Carlo simulations. The sim- 
ulations were performed on a lattice of size 10 and averaged 
2 x 10 7 times with J = AD. 



$2 is a scaling function of the variables rrii/mj and the 
parameters AV^, too and A. In the limit when too — > 
and A — > 00, we expect the AV dependence of $2 to be 
$ 2 (AV) ^ x) . We then conclude that 



7(n) = 



/2rf + 2 
V d+2 



d+2 



n(n — 1) 



0(e 2 ). (13) 



The first term coincides with ^koim [see Eq. QJ]. The 
presence of the second term leads to breakdown of self 
similarity. This is due to the effective anticorrelation be- 
tween the particles. The self similarity conjecture for 
MM is equivalent to a renormalized mean field theory in 
which only the coupling renormalization is taken into ac- 
count. Therefore, it cannot take into account correlations 
between particles. 

Though, Eq. I|13[) is an e-expansion, it is still possi- 
ble to confirm that 7(n) 7^ "/koim(n) in d < 2 by com- 
puting 7(2) exactly. From the definition of 7(2), it fol- 
lows that (R fil R fJ _ 2 ) = (^i^2) 7(2 ^ 2_1 V'(Mi/A 1 2), where 
ip(x) is an unknown scaling function with the property 
ij){x) — i(j(1/x). We need to know the /ii,/Z2 - * be- 
havior of (i? Ml i? M2 ). Averaging Eq. (@J with respect to 
noise and setting dt(Rp) — in the large time limit, we 
find that (R^R^) = j/(Xmo) ~ J/i/A for fi <C tUq . 
Comparing this result with the above scaling form, we 
find that 7(2) = 3 exactly, which coincides with Eq. i|13|) 
with terms of 0(e 2 ) and higher order set to zero. 

Note that 7(2) does not depend on the dimension d. 
The result 7(2) = 3 is a counterpart of the 4/5-th law 
of Navier-Stokes turbulence and is due to conservation 
of mass flux. Recall that the 4/5 law states that, in the 
inertial range the third order structure point function 
of velocity field scales as the first power of separation. 
While the Kolmogorov theory respects 4/5-law in Navier- 
Stokes turbulence, it violates 7(2) = 3 law in MM. 



The terms 0(e 2 ) and higher in Eq. (|13fl are either very 
small or equal to zero when d = 1. We showed their 
absence for n = 1, 2 by means of an exact computation. 
For n = 3, 4 we verify their smallness in one dimensions 
using Monte Carlo simulations (see Fig.^J). 

In two dimensions, logarithmic corrections to the Kol- 
mogorov scaling are expected. These can be calculated 
exactly using the RG method. We present the final re- 
sults here: 



C„(m)> 



J n / 2 [ln(m)]"- n2 /2 



1+0 



1 



In (to 



d=2. (14) 



The logarithmic corrections vanish for n = 2. This is 



consistent with the exact result for C2 (to) 



,-3 



, which 



is valid in all dimensions. 

It is straightforward to extend the above results to CM. 
Kolmogorov theory for the charge model is developed 
using the flux of charge squared J c as a self-similarity 
parameter. This leads to Kolmogorov spectrum for CM, 



CM 



Ikol 



(iTf? ) n - Tne stochastic equation for CM is 



the same as Eq. ||3J modulo a redefinition of the source 
term: j /mo — > 2J C (1 — cosh(/xmo))/mg). Therefore the 
RG analysis presented for MM goes through with just 
one modification: the scale of diffusive fluctuations in 
the charge model is Ljj — ( JcZ) -1 /! -2 )*^ 2 " for /j, — * 0. 
The result of the analysis is 



l CM (n) 



lkohn(n) 



d + 2 



n{n- 1) + 0(e 2 ). (15) 



In complete analogy with MM case, 0(e 2 ) terms in the 
above answer vanish for n — 1 and n — 2. Kolmogorov 
theory correctly predicts the scaling of average charge 
distribution in CM, but fails to predict the scaling of 
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multi point probability distributions of charge. Rela- 
tion 7 CM (2) = 4 reflects conservation of flux of charge 
squared and is satisfied in all dimensions. 

To summarize, we considered two models of aggrega- 
tion with input which had features qualitatively similar 
to turbulence. For these models, we were able to com- 
pute the statistics of the local mass (charge) distribu- 
tion and compare the results against the predictions from 
Kolmogorov theory. It turned out that Kolmogorov self- 
similarity conjecture is equivalent to a renormalized mean 
field theory approximation which correctly accounts for 
effects of reaction rate renormalization in both MM and 
CM. As a result, the answers for average mass (charge) 
distribution derived from Kolmogorov theory turned out 
to be correct. However, Kolmogorov theory failed to 
predict multi point distributions of mass (charge) cor- 
rectly: small scale correlations between diffusing, coagu- 
lating particles are relevant for the multi point statistics, 
but cannot be accounted for by coupling constant renor- 
malization alone. 
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